Excitation and relaxation in atom— cluster collisions 
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Electronic and vibrational degrees of freedom in atom-cluster collisions are treated simultaneously 
and self-consistently by combining time-dependent density functional theory with classical molecu- 
lar dynamics. The gradual change of the excitation mechanisms (electronic and vibrational) as well 
as the related relaxation phenomena (phase transitions and fragmentation) are studied in a common 
framework as a function of the impact energy (eV. . . MeV). Cluster "transparency" characterized 
by practically undisturbed atom-cluster penetration is predicted to be an important reaction mech- 
anism within a particular window of impact energies. 

PACS numbers: 31.70.Hq, 31.15.Ew, 34.50.Bw, 36.40.Ei 
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Collisions with atomic clusters represent a relatively 
new branch of collision physics as compared to the well 
established fields of ion-atom collisions M and ion-solid 
interaction M . The study of cluster collisions is of par- 
ticular interest and importance because it offers the pos- 
sibility to tackle bridge-building questions (like the con- 
tinuous transition from individual excitations in the ele- 
mentary ion-atom collision to the macroscopic stopping 
power in solids) as well as fundamental problems (like 
phase transitions in finite systems). It is also an ex- 
tremely challenging and complicated field as the com- 
prehensive understanding of these collisions still requires 
the development of basically new techniques for both, 
large-scale (multi-parametric) experiment and many- 
body (quantum-mechanical) theory. In this request, the 
present situation resembles very much that of nuclear 
physics at the beginning of the 80's j|. 

Experimentally, large progress has been made, mean- 
while, in the investigation of adiabatic cluster collisions 
where the reaction mechanism is determined by vibra- 
tional excitations only. Typical examples are the study 
of the vibrational energy transfer B , the fusion between 
clusters ||, the formation of endohedral complexes ||, 
and the collision induced dissociation (CID) M. There is 
also a lasting interest to study non-adiabatic cluster col- 
lisions where electronic transitions occur. Experiments in 
this field concern the measurements of the charge trans- 
fer flq-fH]] , ionization and electronic excitation pT| , as well 
as the selective observation of vibrational and electronic 
excitations p"2| . 

Theoretically, adiabatic cluster collisions can be well 
described by quantum molecular dynamics (QMD) or 
molecular dynamics (MD) [|13| — |lT|| . Also, isomeric |lq | 
and solid-liquid phase transitions |19 2fJ in clusters as 
well as the fission process of clusters [|2l[ have been stud- 
ied with MD or QMD where, basically, electronic excita- 
tions are not considered. On the other hand, electronic 
transitions in non-adiabatic cluster collisions have been 



treated with classical [E2 231, semi-classical 24 , or one- 
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electron quantum mechanical |25| , |2q | approaches where 
the atomic structure, and thus the vibrational degrees of 
freedom, are not taken into account. Recently, a gen- 
eral theory has been developed J27| which is able to de- 
scribe simultaneously adiabatic and non-adiabatic colli- 
sions and, in particular, also the still completely unknown 
transition regime where both — electronic and vibra- 
tional — excitations occur. This so-called non-adiabatic 
quantum molecular dynamics (NA-QMD) [G7[ can also 
be used to study, for the first time, phase transitions 
in and fragmentation of clusters induced by electron- 
vibration coupling. 

In this work, different excitation mechanisms (elec- 
tronic and vibrational) as well as related relaxation 
phenomena (phase transitions and fragmentation) in 
cluster collisions are studied in the microscopic frame- 
work of NA-QMD. This theory treats electronic and 
vibrational degrees of freedom simultaneously and self- 
consistently in atomic many-body systems by combining 
time-dependent density functional theory |28] with clas- 
sical MD. The key point, in order to derive practicable 
equations of motion, is the split of the total electronic 
density into the adiabatic and a remaining part, treating 
afterwards the exchange-correlation terms with the adia- 
batic density only [E7J . This approximation has been ap- 
plied and successfully tested against experimental data, 
so far, for the interpretation of fragment correlations in 
CID [G9| and the calculation of absolute cross sections for 
charge transfer pfl] in cluster reactions. Here the univer- 
sal NA-QMD approach is used to investigate the gradual 
change of the excitation and relaxation mechanisms as a 
function of the impact energy in a wide range for two, 
basically different collision systems, Nag + + Na (with at- 
tractive adiabatic forces between projectile and target) 
and Na9 + + He (with dominating repulsive forces and a 
"magic" initial electronic configuration). 

In Fig. 1, the total kinetic energy loss (TKEL) of 



the collision, AE=E cm — E cm (t^> + oc), with E cm the 
impact energy and E cra (t— ► + oo) the final kinetic en- 
ergy of the relative motion between cluster-projectile and 
atomic target in the center-of-mass system is shown cal- 
culated for a fixed collision geometry (with impact pa- 
rameter 6=0) but in a wide range of impact energies 
£ cm =0.2eV...lMeV. 
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FIG. 1. Total kinetic energy loss of relative motion AE 
in central Nag + + Na collisions for a fixed collision geometry 
(see insert) as a function of the center-of-mass impact energy 
E cm calculated with NA-QMD (solid line) and QMD (dotted 
line) . Vibrational and electronic contributions to AE are dis- 
tinguished by grey and light-grey shaded areas, respectively. 
The dot-dashed line corresponds to AE — E cln . 

For "real" scattering events (see below) , AE describes 
the change of the internal energy of the system gen- 
erally connected with excitations in the cases studied 
here. In order to decide between electronic and vibra- 
tional contributions to AE the NA-QMD results are 
compared to those obtained from adiabatic QMD |27| ] 
calculations performed with the same collision geome- 
try. At impact energies below E cm ~ 200 eV, QMD 
and NA-QMD results are clearly seen to be identical 
and, thus, only vibrational excitations occur (adiabatic 
regime) . Above E cm ~ 5 keV electronic transitions dom- 
inate (non-adiabatic regime). In the intermediate range 
of E cul both mechanisms compete (transition regime). 
Maximal TKEL, namely AE=E cm , is practically real- 
ized at all impact energies below E cm ~ 10 eV and is con- 
nected with the formation of relatively long-living but 
in general unstable intermediate compounds Naio + [ pjj l- 
The "real" scattering events with incomplete energy loss 
appear above E cva ~ 10 eV. 

In Fig. 2, the time dependence of AE(t)=E cm — E cm (t), 
with E cm (t) the actual kinetic energy of the rel- 
ative motion, and the displacement of the clus- 
ter atoms d(t) are shown for characteristic impact 
energies. The quantity AE(t) gives insight into 



the collision and excitation dynamics (forces, inter- 
action time, TKEL). The displacement, defined as 
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d(t):=j:^ =1 [Jl A (t)-R A (0)Y 
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(NR) with R A (0) 



the equilibrium positions of the TV atoms in the ground 
state configuration and R the cluster radius, character- 
izes quantitatively the relaxation, i. e., if d is well below 
one a "solid" configuration is described, if d lies in the 
vicinity of one a "liquid" state is realized, and if d goes 
to infinity fragmentation occurs. 
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FIG. 2. Difference of the impact energy and the actual ki- 
netic energy of relative motion AE(t) (left column) and dis- 
placement of the cluster atoms d(t) (right) as a function of 
time t for four impact energies E C m=0.2, 2, 20, 200 keV and 
the same collision geometry as used in Fig. 1, calculated with 
NA-QMD (dark solid lines). In the left column, the corre- 
sponding adiabatic QMD calculations (dotted lines) are also 
shown for comparison. Note the different time scales for the 
excitation processes (femtoseconds) and the relaxation phe- 
nomena (picoseconds). 

Excitation and relaxation mechanisms and, in partic- 
ular, the related time scales are basically different in all 
cases. In the pure adiabatic regime (i? cm =200eV), the 
vibrational excitation of the cluster (A2?«1.4eV) is im- 
mediately connected with a fragmentation process which, 



in contrast to statistical evaporation, starts to proceed 
already during the interaction time of ~40fs (see left 
part of Fig. 2) as a consequence of sufficient momentum 
transfer between projectile and target atoms. Such non- 
statistical fragmentation has been experimentally veri- 
fied only very recently p2| (called "impulsive" fragmenta- 
tion mechanism) . On the contrary, in the high-energetic 
non-adiabatic region £' cm =200keV electronic excitation 
(~ 1 fs) and fragmentation (~ 300 fs) processes are sep- 
arated by two orders of magnitude. Obviously, electron- 
vibration coupling needs time to become effective and 
to induce dissociation ("electronic" fragmentation). At 
an impact energy of E cm =2 keV the cluster remains sta- 
ble after the collision with a total excitation energy of 
AEm0.7 eV originally stored predominantly in electronic 
excitation (cf. adiabatic and non-adiabatic contributions 
to AE in the left part of Fig. 2). Electron-vibration 
coupling leads, however, to the excitation of collective 
surface vibrations indicated by the regular oscillating 
behaviour of d(t) well developed after about 1 ps. At 
E" cm =20keV the first step of the relaxation phase is very 
similar to that observed at _E cm =2keV. But, in this case 
electron-vibration coupling induces clearly a phase tran- 
sition from "solid" to "liquid" in a second relaxation step 
coming about 10 ps after the collision. Owing to the large 
TKEL (AZ?«2.7eV) one should expect statistical evap- 
oration as the final decay channel. 
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FIG. 3. Total kinetic energy loss AE in Na g + + Na col- 
lisions (i? cm =100keV) for randomly chosen cluster orienta- 
tions as a function of the impact parameter b calculated with 
NA-QMD (dots) as well as the resulting mean value (solid 
line). The prediction of the Lindhard model f32l is also shown 
(dashed) . 

The excitation and relaxation mechanisms, summa- 
rized in Figs. 1 and 2, appear in a well defined order 
as a function of the impact energy due to the deliber- 
ately fixed collision geometry. For a given impact energy 
and, in particular, for comparison with experiment the 
impact parameter dependence of measurable quantities 
averaged over different orientations of the cluster with 
respect to the beam axis is of interest. 
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FIG. 4. Upper part: Total kinetic energy loss AE in cen- 
tral Na<) + + He collisions as a function of the center-of-mass 
impact energy E cm calculated with NA-QMD using randomly 
chosen orientations (dots) . The lines correspond to a complete 
energy loss AE — E cul (dashed) and the maximal energy loss 
for a binary collision of the He atom with one cluster atom 
(dotted) |3]IJ33j. Lower part: Total fragmentation cross sec- 
tions as a function of E cnl calculated with NA-QMD (solid 
line). Available experimental data (circles) are from (K3J. 

For the TKEL this dependence is shown in Fig. 3 for an 
impact energy in the non-adiabatic regime. The results 
are compared to the well-known Lindhard model p2j of 
the electronic stopping power in solids. This phcnomcno- 
logical model describes surprisingly well the mean value 
of AE for impact parameters smaller than the cluster 
radius, b ~ 8 a.u. There are, however, two character- 
istic and measurable peculiarities of AE(b) in the mi- 
croscopic calculations. First, for small impact parame- 
ters b ~ 6 a.u., there are a few events with a very large 
TKEL. They result from binary atom collisions between 
the Na-target and a cluster atom, leading to very fast 
fragment monomers. In a macroscopic language, these 
direct "knock-out" collisions correspond to the nuclear 
stopping power H. Second, for impact parameters larger 
than the cluster radius, b ~ 8 a.u., the calculated (in that 
case purely electronic) TKEL are not small as compared 
to the dissociation energy of the cluster. Consequently, 
"electronic" fragmentation (cf. Fig. 2) may essentially 
contribute to the total fragmentation cross section which 
can be expected to be several times larger than the geo- 
metrical one. 

The most interesting and unexpected phenomenon of 



the studies performed so far is the occurrence of clus- 
ter "transparency" in the closed shell system Nag + + He. 
This mechanism is characterized by a practically undis- 
turbed atom-cluster penetration in a certain, small range 
of E cm within the transition regime. In Fig. 4, the TKEL 
as a function of E CTa calculated for random orientations 
with impact parameter b=0 is shown. It can be clearly 
seen that in the window between E cm ~ 10 eV and ~ 
100 eV the variance of AE becomes markedly large and 
there are a lot of events with an extremely small energy 
loss. Obviously, for specific collision geometries in this 
range of E cm the relative velocity between projectile and 
target becomes too large to excite vibrational degrees of 
freedom and is, on the other hand, still too small to in- 
duce electronic transitions. As a consequence, the cluster 
appears to be transparent even in central collisions. 

One possibility to study the change in the excitation 
mechanisms and in particular the occurrence of trans- 
parency is to consider the total fragmentation cross sec- 
tion as a function of E cm (lower part of Fig. 4). The the- 
oretical cross section <7f rag has been calculated by taking 
into account an initial temperature of the clusters accord- 
ing to recent experiments J33[ |. An excellent agreement 
between theory and experiment is found in the adiabatic 
regime where experimental data are available at present, 
i.e. E cm ~ 4eV. It is just the transparency effect which 
leads to a pronounced plateau in crf rag between 10 eV and 
100 eV before electronic excitations induce a further in- 
crease of OjErag(2? cm ). 

In summary, we have shown how different excitation 
mechanisms compete and various relaxation phenomena 
occur in atom-cluster collisions. First relations and basic 
differences to macroscopic systems have been discussed. 
Cluster transparency has been predicted to be an impor- 
tant reaction mechanism in Nag + + He. A systematic 
investigation of different collision systems (in particular 
the size dependence of transparency |34| , its eventual re- 
lation to the channeling effects in solids |p5| ; the tem- 
perature dependence, etc.) will be the scope of future 
studies. 
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